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The thermal reactivity and kinetics of five coal chars, a biomass char, and two coal/biomass char blends in 
an oxy-fuel combustion atmosphere (30%02-70%CC>2) were studied using the non-isothermal ther- 
mogravimetric method at three heating rates. Fuel chars were obtained by devolatilization in an 
entrained flow reactor at 1273 K under N 2 and C0 2 atmospheres. Three nth-order representative gas 
-solid models - the volumetric model (VM), the grain model (GM) and the random pore model 
(RPM) - were employed to describe the reactive behaviour of the chars. The RPM model was found to be 
the best for describing the reactivity of the high rank coal chars, while VM was the model that best 
described the reactivity of the bituminous coal chars, the biomass char and the coal-biomass blend char. 
The kinetic parameters of the chars obtained in N2 and CO2 in an oxy-fuel combustion atmosphere with 
30% of oxygen were compared, but no relevant differences were observed. The behaviour of the blend of 
the bituminous coal (90%wt.) and the biomass (10%wt.) chars resembled that of the individual coal 
concealing the effect of the biomass. Likewise, no interaction was detected between the high rank coal 
and the biomass chars during oxy-fuel combustion of the blend. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Whereas power generation from CCVneutral renewable fuels is 
being promoted, coal utilization is expected to continue in the 
future, as the reserves of coal are abundant and its cost is relatively 
low [1 ]. However, the use of coal in power plants generates a large 
amount of CO2, which is the chief contributor to global climate 
change. One of the promising technologies being developed to 
facilitate carbon capture and storage (CCS) from fossil-fuel-fired 
power plants is the oxy-fuel combustion. The oxy-coal combus¬ 
tion together with carbon sequestration is a strategy that has 
received considerable attention as a viable, cost-effective tech¬ 
nology for power production and that can be used as a retrofit on 
existing coal-fired power plants. This technology involves the 
combustion of pulverized coal in a mixture of oxygen and recycled 
flue gas (RFG), which consists mainly of CO2 and H2O [2], In this 
process, the net volume of flue gas is reduced and a C0 2 volume 
fraction in the flue gas of 95% or higher is provided, which after 
purification can be directly stored in a supercritical state by means 
of compression. 
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The oxy-coal combustion is expected to be different from the 
coal combustion under an O2/N2 atmosphere as the CO2 has higher 
specific molar heat capacity and varying radiation properties 
compared to those of N2. Previous research studies have focused on 
the fact that the use of CO2 instead of N2 may cause a reduction in 
the propagation speed, flame stability and gas temperature profile 
or lead to an increase in the unburned carbon content. Thus, Liu 
et al. [3] concluded that during the oxy-fuel combustion process, 
these problems can be overcome by increasing the oxygen 
concentration in the oxidizer (up to approximately 30%) to match 
the combustion performance achieved in air. However, the reaction 
models and kinetics to describe coal combustion in O2/CO2 atmo¬ 
sphere have been scarcely studied. The knowledge of the effect of 
the CO2 atmosphere on the reactivity of coal is needed for evalu¬ 
ating the kinetic parameters required for computational fluid 
dynamics (CFD) calculations [4], which will be useful for designing 
and modelling oxy-fuel combustion at industrial scale. There are 
a few recent studies on the calculation of the kinetic parameters of 
the combustion of coal in an oxy-fuel atmosphere. Niu et al. [5,6] 
carried out the coal combustion in O2/CO2 mixtures on a thermog- 
ravimetric analyzer and calculated the kinetic parameters which 
describe the combustion process by linearization using the Coats- 
Redfern method. Wang et al. [7] compared the combustion of coals 
and chars in oxy-fuel atmosphere applying the non-isothermal 
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thermogravimetric analysis, but these authors claimed that insuf¬ 
ficient information exits on the effect of char preparation atmo¬ 
sphere on the char reactivity. 

Thermogravimetric analysis (TGA) is a common technique used 
to investigate and compare thermal events and determine kinetic 
parameters during the combustion, pyrolysis and gasification of 
solid raw materials, such as coal, wood, etc. [8—14], Moreover, 
quantitative methods can be applied to TGA curves to obtain kinetic 
parameters of the thermal events. Miura and Silveston [15] showed 
the validity of the temperature-programmed reaction (TPR) tech¬ 
nique for the analysis of noncatalytic gas—solid reactions. This 
technique appears to provide more kinetic information than what 
is obtainable from the same number of experiments performed at 
constant temperature. Kasaoka et al. [16] also stated that in an 
isothermal experiment, a tedious repetition of experimental runs is 
required to determine the kinetic parameters of the Arrhenius 
equation. 

On the other hand, with the EU announcing that it intends to 
supply 20% of its overall energy needs from renewable sources by 
2020, interest in biomass as a renewable source is growing [17], 
This source of energy is considered carbon neutral because the 
carbon dioxide released during biomass utilization is recycled as an 
integral part of the carbon cycle. Co-combustion of biomass and 
coal has generated widespread interest because of the reduced 
emissions of gases such as CO2, SO2 and NO x compared to those 
emitted by the combustion of coal [18], The cofiring of biomass 
with coal is advantageous since biomass has higher volatile matter 
content and lower devolatilization temperature, which can aid the 
ignition and combustion characteristics of the blend [19], Lai et al. 
[20] compared the combustion behaviour of lignocellulosic mate¬ 
rials in CO2/O2 and N2/O2 atmospheres and calculated the kinetic 
parameters considering three parallel reactions during the 
combustion. A few studies have been published on the oxy-fuel 
combustion of coal and biomass blends [21—23], but reactivity 
and kinetic studies of biomass and coal/biomass blends under oxy- 
fuel combustion needs to be carried out. 

The aim of the present work was to study the oxy-fuel 
combustion (30%02-70%C02) reactivity and kinetic behaviour of 
five coal chars and a biomass char obtained under 100%N2 and 100% 
CO2 atmospheres, as well as the coal/biomass char blends. For this 
purpose, the temperature-programmed reaction (TPR) technique at 
three different heating rates was used. Three mathematical 
models - the volumetric model (VM), the grain model (GM) and 
the random pore model (RPM) — were used to determine the 
kinetic parameters which best represent the oxy-fuel combustion 
characteristics of the coal char, biomass char and coal/biomass char 
blends under an oxygen-carbon dioxide atmosphere. 

2. Experimental 

2.1. Fuel samples 

Five coals were used in this work: an anthracite (AC), a semi¬ 
anthracite (HVN) and three high-volatile bituminous coals (DAB, 
M6N and NZ). A type of biomass, torrefacted pine sawdust (TPIN), 
was also employed. The coal and biomass samples were ground, 
sieved and the resulting 75-150 pm size fraction was used for the 
devolatilization tests. The results of the proximate and ultimate 
analyses and high heating values of the coal and biomass samples 
are shown in Table 1. 

The torrefaction of the biomass was carried out to improve its 
properties for pulverised systems [24], To obtain the TPIN sample, 
a raw pine sawdust sample was ground and sieved to obtain 
a particle size fraction of 75—212 pm. Then, the torrefaction of this 
sample was performed using a vertical quartz reactor with 


Table 1 

Proximate and ultimate analyses and high heating value of the samples. 


Sample 

AC 

HVN 

DAB 

M6N 

NZ 

PIN 

TPIN 

Origin 

Spain 

Spain 

China 

Mexico 

1 New Zealand 

Spain 

Spain 

Rank 

an 

sa 

hvb 

hvb 

hvb 



Proximate analysis* 







Moisture (wt.%) 

2.3 


2.9 

1.8 

11.5 

6.8 


Ash (wt.%, db) 

14.2 

10.7 

10.9 

30.2 

2.9 

3.8 

4.2 

V.M. (wt.%, db) 

3.6 

9.2 

28.8 

30.6 

47.9 

79.8 

75.5 

F.C. (wt.%, db) b 

82.2 

80.1 

60.3 

39.2 

49.2 

16.4 

20.3 

Ultimate analysis 

(wt.%, daf) a 






C 

94.7 

91.7 

81.9 

76.7 

67.5 

45.9 

51.2 

H 

1.6 

3.5 

5.0 

6.2 

5.3 

6.1 

5.7 

N 

1.0 

1.9 

1.1 

1.6 

1.0 

0.7 

0.9 

S 

0.7 

1.6 

1.2 

1.2 

0.2 

0.0 

0.0 

O b 

2.0 

1.3 

10.8 

14.3 

26.0 

47.3 

42.2 

HHV (MJ/kg, db) 

29.2 

31.8 

28.8 

23.1 

27.9 

18.9 

20.2 


an: anthracite; sa: semi-anthracite; hvb: high-volatile bituminous coal, 
db: dry basis; daf: dry and ash-free bases. 

a The proximate analysis was conducted in a LECO TGA-601, and the ultimate 
analysis in a LECO CHNS-932. 
b Calculated by difference. 

a diameter of 20 mm and a length of 50 mm, where 10—15 g of 
biomass was kept at a temperature of 513 K for 1 h under a nitrogen 
flow rate of 50 NmL min -1 . 

2.2. Char preparation 

The chars were prepared by devolatilizing the raw fuels in an 
electrically heated entrained flow reactor (40 mm internal diam¬ 
eter, 1400 mm length) in streams of 100%N2 and 100%CO2 
(4.79 NL min -1 ). The experimental device has been described 
elsewhere [25,26], The devolatilization experiments were carried 
out at a reactor temperature of 1273 K and a particle residence time 
of 2.5 s. After the experiments, the chars were cooled down under 
a flow of nitrogen to room temperature and a water-cooled col¬ 
lecting probe was inserted into the reaction chamber from below to 
collect the char samples. The obtained coal and biomass char 
samples were used for the oxy-fuel combustion reactivity experi¬ 
ments. Moreover, two blends, AC + TPIN and DAB + TPIN, 
composed of 10wt% of biomass char TPIN-CO2 and 90wt% of coal 
char AC-CO2 or DAB-C0 2 , respectively, were prepared. The 
components of the mixtures were blended in the adequate 
proportions and manually homogenized. 

2.3. Oxy-fuel combustion reactivity tests of the chars 

The reactivity tests were conducted in a thermobalance 
(Setaram TAG24) at atmospheric pressure. Approximately 5 mg of 
sample was placed in a crucible of 2 mm height with a circular base 
of 5 mm diameter. A small amount of sample and slow heating rates 
were used to avoid heat transfer limitations and to minimize mass 
transfer effects. A thermocouple was located close to the platinum 
basket to monitor the temperature and to close the oven control 
loop. In this work, all the experiments were performed under non- 
isothermal conditions from room temperature to 1273 K at three 
different heating rates: 2,3 and 5 K min -1 . For the NZ samples, the 
heating rates used were: 1, 2 and 3 K min -1 . The total flow rate of 
the gas introduced into the thermobalance during the oxy-fuel 
combustion experiments was 50 NmL min -1 , the gas consisting of 
30%C>2 and 70%C02. The char conversion, X, and the reaction rate, 
dX/dt, were calculated. The TG reactivity tests were firstly per¬ 
formed for the coal and biomass chars separately, and then for the 
coal/biomass char blends. 
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3. Kinetic models 

A general kinetic expression for the overall reaction rate in gas- 
solid reactions is written as follows [27]: 

dX/dt = k(P g , T)f(X) (1) 

where k is the apparent combustion reaction rate, which includes 
the effect of temperature (T) and the effect of the reactive gas 
partial pressure ( P g ), and J(X) describes the changes in the physical 
or chemical properties of the sample during reaction. X represents 
the loss in mass fraction or mass conversion ratio, which can be 
calculated by the following relationship: 

X = (m 0 - m t )/(m 0 - m f ) (2) 

where mo is the initial mass of the char sample, m t the mass of the 
char sample at time t and mj the final mass of the char sample. 

Assuming that the partial pressure of the reactive gas remains 
constant during the process, the apparent combustion reaction rate 
will be dependent on the temperature and can be expressed using 
the Arrhenius equation, as follows: 

k = koe- £ / Rr (3) 

where ko and E are the pre-exponential factor and activation 
energy, respectively. 

In this work, three nth-order models were applied to describe 
the reactivity of the chars studied: the volumetric model (VM), the 
grain model (GM) and the random pore model (RPM). These 
models give different formulations of the term J[X), with X repre¬ 
senting the degree of char conversion on a dry ash-free basis. 

The VM assumes a homogeneous reaction throughout the 
particle and a linearly decreasing reaction surface area with 
conversion [28], The overall reaction rate is expressed by: 

dX/dt = fcvM(l-X) (4) 

The GM or shrinking core model, proposed by Szekely and Evans 
[29], assumes that a porous particle consists of an assembly of 
uniform nonporous grains and that the reaction takes place on the 
surface of these grains. The space between the grains constitutes 
the porous network. The shrinking core behaviour applies to each 
of these grains during the reaction. In the regime of chemical 
kinetic control and, assuming the grains have a spherical shape, the 
overall reaction rate can be expressed as: 

dX/dt = k GM (l -X) 2 / 3 (5) 

This model predicts a monotonically decreasing reaction rate 
and surface area because the surface area of each grain is receding 
during the reaction. 

The RPM model considers the overlapping of pore surfaces, 
which reduces the area available for reaction [30]. The basic 
equation for this model is: 

dX/dt = W1 — X)[l - i/4n(l -X)] 1 / 2 (6) 

This model is able to predict a maximum value of reactivity as 
the reaction proceeds, as it considers the competing effects of pore 
growth during the initial stages of combustion and the destruction 
of the pores due to the coalescence of neighbouring pores during 
the reaction. The RPM model employs two parameters, the 
apparent combustion reaction rate constant, Icrpm, and \p, which is 
related to the pore structure of the unreacted sample, although 
Miura and Silveston [15] suggested that the i/' parameter should be 
considered as an adjustable parameter. 


The non-isothermal thermogravimetric method or temperature- 
programmed reaction (TPR) technique involves heating the samples 
at a constant rate, a. The temperature, T, is related to time, t, by: 

T=T 0 + at (7) 

where T 0 is the temperature at which heating is started, which can 
be set equal to 0 provided that To is low enough for the reaction rate 
to be practically zero when heating is initiated. 

By means of Eq. (7), Eq. (4) can be integrated to give: 

X= 1-exp (-^p(ti)) (8) 

where 

P(u) = S — ~ J u (9) 


u = E/RT (10) 

From the literature, several proposed approximations for p(u) 
can be found. In this study the one employed has been [15,31,32]: 

p(u) = e u /u 2 (11) 

This approximation is valid for u > 10, which is fulfilled by these 
fuels when burned in the studied oxy-fuel atmosphere. Eq. (8) can 
then be written as: 

X = l-exp[-^] (12) 

Similarly, Eqs. (5) and (6) can be integrated with the above 
approximation, to give Eqs. (13) and (14) respectively: 

(,3) 

x - i - exp [-S t " ei '( ,+ j(S) t " i!if )] (i4) 

Eqs. (12)—(14) can be utilized to determine ko and E values from 
X experimental data sets by employing nonlinear least-squares 
fitting methods. According to Miura and Silveston [15], the deter¬ 
mination of the kinetic parameters from a single TPR run may lead 
to unreliable rate parameters and, furthermore, the fitting of data 
by a model may not validate the model if just one TPR run is used. 
These authors claimed that at least three TPR runs at different 
heating rates are required to estimate reliable parameters and 
accurate activation energies. 

In this study the kinetic parameters were determined from 
three TPR runs, each one performed at a different heating rate. The 
nonlinear least-squares method was employed to fit the experi¬ 
mental data of 1-X vs. temperature, T, to the three models by 
means of Eqs. (12)—(14), and to estimate the ko and E values which 
minimize the objective function, OF: 

0F = E [(! X >exp,i - (1 -*)calcj 2 (15) 

where (1 -X) ex p,i is the experimental point corresponding to the ith 
temperature, T„ (1 —X) ca i Ci ,- is the value calculated at Tj, and N is the 
number of data points. 
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According to Varhegyi [33], the method of least-squares can 
work better in non-isothermal kinetics than other methods 
requiring simpler computer programming or less computational 
time, since this technique aims directly at the description of the 
experimental data in a wide range of experimental conditions. 

To assess the quality of the fit and verify the capacity of the 
kinetic models to describe the degree of char conversion, X (or 
1—X), the deviation (DEV) between the experimental and calcu¬ 
lated curves was obtained using the following expression: 

DEV(1 -X)(%) = 100 ((1 -X) expj 

1=1,N 

- (1 -^)caic,i) 2 /w] 1/2 /max(l -X) exp (16) 

where (1-X) ex p,i and (1-X) ca ic,i represent the experimental and 
calculated data of 1-X, N is the number of data points, and 
max(l-X) exp is the highest absolute value of the experimental 
curve. The best fitting kinetic parameters and model were chosen 
from the highest R 2 value obtained from the results which proved 
to be statistically significant, together with the lowest value of DEV 
(1-X). 

4. Results and discussion 

4.1. Kinetic parameters 

Table 2 shows the kinetic parameters (E, k 0 and f) determined 
from the data obtained at three heating rates (1,2 and 3 K min -1 for 
the NZ coal char and 2, 3 and 5 K min -1 for the rest of samples) 
together with the coefficients of determination, R 2 , for each model 
and sample. R 2 shows the variation in the dependent variable, 1 -X, 
which is explained by the model. Table 2 also presents the statis¬ 
tically significant model fittings. 

The RPM model fits the experimental data better than the other 
two models for coal chars AC-N2 (R 2 = 0.999), AC-CO2 (R 2 — 0.999), 
HVN-N 2 (R 2 = 0.998) and HVN-CO2 (R 2 = 0.998), since it displayed 
a significant fit and presented the highest R 2 value (Table 2). In the 
case of the DAB, M6N and NZ char samples, as well as the 
DAB + TPIN blend, the i/' value obtained from the RPM model is 
a meaningless negative value, which indicates that the best fit to 
the experimental data is obtained with a tp value of zero. For the 
TPIN samples, the \p value obtained from the RPM model is close to 
zero and the fits were not significant. When the i/' value is close to 
zero, the RPM model predicts a nearly constant decrease in reac¬ 
tivity with conversion, as does the VM model. Therefore, it can be 


concluded that the model which best describes the reactivity of the 
DAB-N 2 (R 2 = 0.998), DAB-CO2 (R 2 = 0.998), M6N-N2 (R 2 = 0.995), 
M6N-CO2 (R 2 = 0.996), NZ-N 2 (R 2 = 0.998), NZ-C0 2 (R 2 = 0.996), 
TPIN-N2 (R 2 = 0.998) and TPIN-CO2 (R 2 = 0.996) char samples, as 
well as the DAB + TPIN (R 2 = 0.999) blend, was the VM model 
(Table 2). On the other hand, the reactivity of the AC + TPIN blend 
has not been satisfactorily described by none of the considered 
models, since the oxy-fuel combustion process took place in two 
stages and it was not then possible to obtain a simple set of kinetic 
parameters for this sample. 

The RPM model predicts a maximum in reaction rate as the 
reaction proceeds, as pore overlapping is considered. Pore shape is 
assumed cylindrical and supposed to grow radially while reaction 
proceeds, instead of keeping initial volume. Initially, the cylinder 
growth causes an increase in total reaction surface, which means 
higher reaction rate. Finally, reaction progress brings about 
a neighbouring pore intersection. Due to pore overlapping, the 
reaction surface area is lower and, consequently, the reaction rate 
decreases [34], Conversely, the VM and CM models cannot describe 
a maximum in reaction rate but predict a constant decrease in the 
reaction rate with conversion. Therefore, from the results of the 
model fits it can be deduced that, in the case of the high rank coals 
used in this work, AC and HVN, the char reaction rate presents 
a maximum value during combustion. However, for the bituminous 
coals (DAB, M6N and NZ), the biomass sample (TPIN) and the coal/ 
biomass blend (DAB + TPIN) the char reaction rate decreases 
linearly during the whole oxy-fuel combustion process under the 
conditions of the present study. 

Figs. 1 and 2 show, for the three heating rates, the experimental 
1-X data and the 1-X curves calculated (Eqs. (12)—(14)) using the 
parameters obtained from the data at the three heating rates for the 
meaningful and statistically significant models (Table 2). As 
example, the curves for the AC, HVN and DAB char samples ob¬ 
tained both in N2 and CO2 are shown in Fig. 1 and the curves for the 
biomass char, TPIN, and the blends, AC + TPIN and DAB + TPIN, are 
shown in Fig. 2. To quantify the errors produced by the kinetic 
models in predicting the values of conversion, the deviation (DEV) 
between the experimental and calculated 1 -X curves was obtained 
using the Eq. (16). The results obtained from the significant models 
for all the char samples and the blend are summarised in Table 3. In 
accordance with the previous results (Table 2), the lowest deviation 
from the calculated values of the conversion rates was obtained 
using the RPM model for all the AC and HVN char samples and the 
VM model for all the DAB, M6N, NZ, TPIN and DAB + TPIN samples 
(Table 3). Saravanan et al. [4] and Gil et al. [35] have also studied the 
reactivity of coal chars during oxy-fuel combustion using the 


Table 2 

Kinetic parameters of the chars during oxy-fuel combustion (30%02-70%C0 2 ) determined with the TPR technique at three heating rates for VM, GM and RPM models. 


Char Volumetric model (VM) Grain model (GM) Random pore model (RPM) 

£ (kj mol -1 ) k 0 (s- J ) R 2 Etkjmol 1 ) Ms -1 ) R 2 Etkjmor 1 ) ko(s~’) R 2 

AC-N 2 
AC-CO 2 
HVN-N 2 
HVN-COj 
dab-n 2 
dab-co 2 

M6N-N 2 
M6N-C0 2 
nz-n 2 
nz-co 2 
tpin-n 2 
tpin-co 2 

DAB + TPIN 


0.977 

0.987 

0.995 


.75E+09 

.10E+06 

.82E+07 

.47E+09 

.73E+08 

.13E+07 

.37E+08 

.60E+05 

.29E+06 

.42E+06 


1.83E+09 
3.54E+07 
4.13E+06 
1.25E+07 
2.95E+04 
5.27E+04 
1.83E+04 
3.48E+04 
2.81 E+02 
2.78E+01 
1.31E+04 
3.72E+04 
5.03E+03 


7.72E+03 

2.03E+07 

7.25E+07 

2.77E+09 

1.57E+09 

2.17E+07 

1.60E+08 

2.68E+05 

5.42E+06 

3.73E+06 


0.999 a 

0.999 a 

0.996 a 

0.997 a 


Statistically significant. 













random pore model. Gil et al. [35] found that RPM model was the 
best for describing the reactivity of coal chars. 

Table 4 summarizes the activation energy values obtained for 
the char samples used in the present study, together with data 
from previous works on oxy-fuel combustion. In the present work, 
using the models of best fit, the activation energy for the coal 
chars under 30%02-70%C02 atmosphere was 116-171 kj mol -1 
(Table 2). In a previous study, Gil et al. [35] had calculated acti¬ 
vation energy values for coal chars under oxy-fuel atmosphere 
with 30%02 of 117—127 kj mol -1 . Niu et al. [5] obtained activation 


energy values for pulverized coals in an oxy-fuel atmosphere (20% 
O2) that ranged between 109 and 248 kj mol -1 . Similarly, Liu [36] 
achieved activation energy values for coal chars in an oxy-fuel 
atmosphere (10% O2) which ranged between 115 and 
147 kj mol -1 . On the other hand, the activation energy for the TPIN 
biomass char under 30%02-70%C02 atmosphere was 118— 
129 kj mol -1 (Table 2). In the case of the DAB + TPIN char 
blend, a value of the activation energy of 136 kj mol -1 was ob¬ 
tained, which is an intermediate value between those of the 
individual components of the blend. 






























Fig. 2. Experimental conversion curves for TPIN-N 2 (a), TPIN-C0 2 (b), AC + TPIN (c) and DAB + TPIN (d) samples during oxy-fuel combustion (30%0 2 —70%C0 2 ) and those calculated 
with three nth-order reaction models (VM, GM and RPM) using parameters determined from heating rates at 2, 3 and 5 K min -1 . 


Table 2 shows that the activation energy values and pre¬ 
exponential factors for the two devolatilization atmospheres are 
very similar. From these data, it can be concluded that no differ¬ 
ences in the kinetic parameters of oxy-fuel combustion under 30% 
0 2 -70%C0 2 atmosphere are produced by using 100%N 2 and 100% 
C0 2 devolatilization atmospheres. 

Wang et al. [7] studied the devolatilization of coal samples in 
100%N 2 and 100%CO 2 by means of TGA and found that the TG and 


Deviation (%) between the experimental and calculated conversion (1-X) data. 


AC-N 2 

AC-C0 2 

hvn-n 2 

hvn-co 2 

dab-n 2 

dab-co 2 

M6N-N 2 

M6N-C0 2 


tpin-n 2 

tpin-co 2 

DAB + TPIN 


DEV (1—X) (%) 

VM GM 

6.39 7.94 

4.36 5.40 

2.80 3.07 

2.30 2.72 

1.53 3.45 

1.63 4.09 

2.89 6.73 

2.52 6.11 

1.63 8.01 

2.67 11.84 

1.82 2.16 

2.30 3.18 

1.08 3.50 


RPM 

1.24 

1.08 

1.57 

1.41 


DTG curves in both atmospheres were almost superposing in the 
low temperature zone (up to ~1073 I<), concluding that the 
pyrolysis reactivity is highly similar in N 2 and C0 2 . At higher 
temperatures, however, in C0 2 atmosphere char-C0 2 gasification 
was observed, which is also in accordance with the results obtained 
by Rathnam et al. [37], At the temperature used in the present study 
for obtaining the chars (1273 K), gasification process could there¬ 
fore occur at some extent, which would affect the char character¬ 
istics. The additional char-C0 2 reaction could then increase the 
reactivity of the coal char obtained under C0 2 atmosphere. 
Nevertheless, the endothermic effect of char-C0 2 reaction could 
also play a negative role in the enhancement of the reactivity 
through reducing the coal particle temperature. Wang et al. [7] 
pointed out the necessity of studies to clarify the circumstances 
in which those opposite effects will preponderate the other. These 
authors concluded that char reactivity (in 21%0 2 /79% C0 2 ) prepared 
in C0 2 is almost equivalent to that prepared in N 2 in a study where 
coal chars were obtained at low temperature (~1073 K) in 
a tubular furnace. Borrego et al. [38] studied the characteristics of 
the chars from various biomass materials generated at a high 
heating rate in a drop tube furnace (DTF) under N 2 and C0 2 
atmospheres and they found that similar morphology, optical 
texture and specific surface area were shown by the biomass chars 
generated under both atmospheres. 

The reaction rate of a char under an oxygen atmosphere is 
usually very high and in the present study the reactivity was 
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Table 4 

Comparison of the activation energy values obtained in the present work and in previous oxy-fuel combustion studies. 


Sample 
Coal chars 
Coal chars 
Pulverized coals 
Coal chars 

Coal/biomass 
char blend 


Atmosphere 
30%O 2 —70%CO 2 
30%O 2 —70%CO 2 
20%O 2 —80%CO 2 
10%O 2 —90%CO 2 
30%O 2 —70%CO 2 
30%O 2 —70%CO 2 


E (kj mol -1 ) 

116- 171 kjmor 

117- 127 kj mor 
109-248 kj mol- 
115-147 kj mol- 

118- 129 kjmor 
136 kj mor 1 


nth-order gas-solid models 
nth-order gas-solid models 
Coats-Redfern method 
Isoconversional model-free methods 
nth-order gas-solid models 
nth-order gas-solid models 


Reference 
Present study 
Gil et al. (2012) [35] 
Niu et al. (2009) [5] 
Liu (2009) [36] 
Present study 
Present study 


measured in a high oxygen content (30%) atmosphere, which may 
have concealed possible differences in reactivity between both 
chars and may explain why no differences were observed 
between the chars obtained under N2 or CO2. Varhegyi et al. [39] 
studied the coal char kinetics in 02/Ar and O2/CO2 atmospheres 
by non-isothermal thermogravimetry and concluded that the 
char reaction rate was proportional to the O2 concentration of the 
ambient gas and was not influenced by the presence of high 
amounts of CO2 since char-C02 reactions have much lower rates 
than char oxidation. Niu et al. [6] confirmed that for O2 concen¬ 
trations of 30—40% in oxy-fuel atmospheres, the reaction process 
is mainly dominated by the diffusion of the O2 to the surface of 
the solid coal and the effect of the carrier gas is not very 
important. 

The relative reactivity of O2 with coal is much greater than CO2 
with coal, which suggests that the char oxidation in the present 
work was predominantly due to oxygen and a possible gasifica¬ 
tion reaction due to CO2 during the reactivity experiments would 
have not existed under these conditions. It was also noted by 
Saravanan et al. [4] in their study about the reactivity of coal chars 
during oxy-fuel combustion under different concentrations of 
C0 2 , who pointed out that the effect of CO2 concentration on the 
char oxidation rate in the presence of oxygen might be only due 
to the bulk resistance offered by CO2 during the char-oxygen 
reaction, although this fact would require more additional 
studies. 

4.2. Comparison of isothermal and non-isothermal reactivity tests 

In a previous study [35], a kinetic analysis of the oxy-fuel 
combustion of the HVN-N2 and HVN-CO2 char samples had been 
carried out by the authors at constant temperature. The results 
obtained using the TPR technique were then compared with those 
obtained from experiments performed at constant temperature. 
From the isothermal oxy-fuel combustion experiments, it was 
concluded that the best model for describing the behaviour of both 
HVN samples was RPM and the parameters estimated using this 
model by both techniques were therefore compared. Following the 
indications from Miura and Silveston [15], the values of koe '^ Rr 
were calculated using the kinetic parameters in Table 2. They were 
plotted on an Arrhenius diagram (Fig. 3) and compared with those 
obtained in the isothermal experiments [34], A good agreement can 
be observed between the k 0 e _£/RT values estimated by both 
methods, indicating that the TPR technique provides reliable 
kinetics parameters when data from the three heating rates are 

4.3. Interactions between the components of the coal/biomass 
blends 

Fig. 4 shows the experimental dX/dt curves for the coal/biomass 
blends (AC + TPIN and DAB + TPIN) under the oxy-fuel combustion 
atmosphere studied, together with the theoretical reaction rate 


curves calculated according to the additive rule from those of the 
individual components. The dX/dt curves for the individual 
components of the blends were also included. The theoretical and 
experimental dX/dt curves of the coal/biomass blends were 
compared to find out whether the components of the blends 
interacted during the oxy-fuel combustion process. As mentioned 
above, the theoretical dX/dt curves of the blends were calculated 
according to the additive rule of blends, i.e.: 

(dX/dt) blend = Xi(dX/dt) coal +x 2 (dX/dt) biomass (17) 

where (dX/dt) coa i and (dX/dt) biomass are the reaction rate of the 
individual fuels, and x\ and X2 are the proportions of coal and 
biomass in the blend, respectively. 

Fig. 4a shows that when the AC high rank coal and the TPIN 
biomass are co-combusted under oxy-fuel atmosphere, two inde¬ 
pendent peaks are observed in the DTG curves of the blend. 
Furthermore, no significant deviations are appreciated between the 
experimental and theoretical dX/dt curves in the case of the 
AC + TPIN char blend at the three heating rates. This suggests that 
no significant interactions occur during char oxy-fuel combustion 
with 30% of oxygen, reflecting the additive behaviour of blend 
components during the combustion process. This means that it 
should be possible to predict the experimental reactivity curve of 
the blend, on the basis of the experimental reactivity curves of each 
individual component and their percentages in the blend. The 
absence of synergetic effects during the char oxy-fuel combustion 



1 1T (1C 1 ) 

Fig. 3. Comparison between the apparent oxy-fuel combustion reaction rates for HVN- 
N 2 and HVN-C0 2 obtained from TPR data (at heating rates of 2, 3 and 5 K min -1 ) and 
data obtained from isothermal experiments [35]. 
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Fig. 4. Comparison between the experimental and theoretical (according to the 
additive rule from those of the individual components) reaction rate curves for the 
AC + TPIN (a) and DAB + TPIN (b) blends during non-isothermal (2,3 and 5 K min -1 ) 
oxy-fuel combustion. 


process, indicates that the oxy-fuel combustion reactions of the 
coal char did not seem to be influenced by the presence of 10wt% of 
biomass char. Each component of the mixture behaved indepen¬ 
dently and did not interact with the other material. 

When the experimental and theoretical dX/d t curves for the 
DAB + TPIN blend are compared (Fig. 4b), it can be seen that dX/dt 
curve of the blend closely resembled that of the coal char sample. 
This indicates that the addition of the TPIN biomass char to the DAB 
coal char, when added in a proportion of 10%, has very little effect 
on the behaviour of the individual coal char. 

5. Conclusions 

The chars from five coals (AC, HVN, DAB, M6N and NZ) and 
a torrefacted biomass (TPIN) obtained in N2 or CO2 atmospheres, as 
well as two coal/biomass blends (AC + TPIN and DAB + TPIN), were 
subjected to oxy-fuel combustion in a thermobalance to study their 


thermal reactivity. The temperature-programmed reaction tech¬ 
nique employed in the analysis of noncatalytic gas—solid reactions 
was applied at three heating rates to estimate the kinetic param¬ 
eters which best described the reactive behaviour of the chars 
during oxy-fuel combustion. The best model for describing the char 
oxy-fuel combustion of the AC and HVN high rank coals was found 
to be the RPM model, whereas the VM model better described the 
char reactivity of the DAB, M6N and NZ bituminous coals, the TPIN 
biomass and the DAB + TPIN blend under oxy-fuel conditions. The 
kinetic parameters showed that no differences in reactivity under 
a typical oxy-fuel combustion atmosphere (30%C>2—70%C02) exist 
between chars obtained under CO2 and N 2 . No significant interac¬ 
tions were detected between the components of the coals-biomass 
blends chars (90-10wt%) with 30% of oxygen. 
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